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We use three Chandra observations, totalHng 134.3 ks, to constrain the dynamical 
motion of NGC 1404 falhng towards the dominant elhptical NGC 1399 through the 
Fornax cluster gas. The surface brightness profile of NGC 1404 shows a sharp edge 
at ~ 8kpc from its center in the direction of NGC 1399, characteristic of jump-like 
, temperature and density discontinuities from ram-pressure stripping of the galaxy gas 

\ caused by its motion through the surrounding ICM. We find the temperature of the 

galaxy gas inside the edge is ~ 2.8 times cooler {kT = 0.55l[J;j]^ keV with abundance 
O ; ^ = O.73l[5:?6^0) than the cluster gas {kT = 1.53lg:J^ keV, ^ = O.42+°;23Z0). We 

■ use the shape of the surface brightness profile across the edge to fit the position of 
the edge, the power law behavior of NGC 1404's density distribution in the leading 

Q_i' direction, and the density discontinuity at the edge. The electron density inside the edge 

Q . (3.9 — 4.3 X 10~^ cm~^ ) depends strongly on the gas abundance; while the density of the 

|H \ ICM (7 — 8 X 10~^ cm~^ ) depends strongly on the assumed geometry (relative distance) 

^ \ between NGC 1404 and NGC 1399. The corresponding pressure jump of 1.7 — 2.1 across 

the leading edge of the galaxy and the cluster free-stream region implies near sonic 

■ motion (Mach number 0.83 — 1.03) for NGC 1404 with velocity 531 — 657 km s~^ relative 

■ to the surrounding cluster gas. The inclination angle of the motion inferred using the 
^ relative radial velocity between NGC 1404 and NGC 1399 as representative of that 

between NGC 1404 and the cluster ICM, is uncomfortably large (> 40°) given the 
sharpness of the surface brightness edge, suggesting either a non-zero impact parameter 
between NGC 1404 and NGC 1399 or that NGC 1399 is also moving radially with 
respect to the cluster ICM. 

Subject headings: galaxies: clusters: general - galaxies: individual (NGC 1404) - galax- 
ies: intergalactic medium - X-rays: galaxies 



1. Introduction 



Our early picture of galaxy clusters, as simple, dynamically relaxed systems, has been changed 
dramatically over the past decade by X-ray images taken with the Einstein Observatory, ROSAT, 
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BBXRT, and ASCA. These images show that galaxy clusters are complex systems with extensive 
substructure and mergers on time sales of ~ 10^ years, (see e.g. Forman et al. 2001, 2002 and ref- 
erences therein) . The excellent angular resolution of the Chandra X-ray Observatory made possible 
the detailed study of these substructures, particularly near the centers of clusters, revealing surpris- 
ing new signatures for their interactions. Notably ram pressure stripping of subcluster gas due to 
its motion through the intracluster medium (ICM) was found to produce prominant X-ray surface 
brightness discontinuities from cold fronts (dense surviving subcluster cores), rather than from the 
shocked gas also expected from such motion (Forman et al. 2001, 2002; Vikhlinin et al. 2001; 
Markevitch et al. 2000; Heinz et al. 2003). Furthermore, high quality X-ray measurements of the 
geometry and gas properties along the interface between the cold front and the ICM were used 
to determine the dynamical motion of the cold front through the cluster (Vikhlinin et al. 2001; 
Mazzotta et al. 2001). 

Galaxies in these dynamically rich cluster environments are subject to both tidal and hydrody- 
namical interactions that significantly affect their evolution (see, Gnedin 2003; Acreman et al. 2003 
and references therein). The Chandra X-ray Observatory has also made possible detailed studies 
of hot gas in member galaxies of nearby clusters in interaction with each other and their gaseous 
environment. Just as in the case of larger cluster substructure, the action of ram pressure acting 
on gas within an individual galaxy due to the galaxy's motion through the ICM produces sharp 
X-ray surface brightness discontinuities along the galaxy's leading edge associated with cold fronts, 
as well as trailing wakes and comet-like debris tails of gas swept from the infalling system (e.g. see 
Forman ct al. 1979, White et al. 1991, Rangarajan et al. 1995 for M86, and Irwin & Sarazin 1996 
for NGC 4472, both in Virgo; and Wang et al. 2004 for C153 in Abell 2125). Thus, in addition to 
elucidating the hydrodynamics of ram pressure stripping itself, high quality X-ray measurements of 
the geometry of these features combined with measurements of the gas temperatures, abundances, 
and densities along the interface between the cold front and the ICM may be used to constrain the 
three dimensional (transverse as well as radial) motion of the galaxy through the cluster (Dosaj 
et al. 2002). 

In this paper we present this analysis for one such galaxy, NGC 1404 (shown in Figure 1), 
that is undergoing ram pressure stripping as it falls inwards towards the Fornax cluster center in 
the direction of the dominant elliptical galaxy NGC 1399. Both NGC 1399 and NGC 1404 are 
well studied in X-rays. Early observations of NGC 1399 using EXOSAT, Einstein IPC, Ginga, and 
BBXRT (see Jones et al. 1997 and references therein) established the presence and mean tem- 
perature of its X-ray emitting halo. More recent measurements by ASCA, ROSAT, XMM-Newton 
and Chandra mapped the temperature and density profiles for the galaxy and the surrounding 
Fornax ICM (Jones ct al. 1997; Buote 2002; Paohllo et al. 2002; O'Sullivan et al. 2003; Scharf 
et al. 2004). NGC 1404 has been observed with ASCA (Loewenstein et al. 1994; Buote & Fabian 
1998), ROSAT PSPC (Jones et al. 1997; Paolillo ct al. 2002; O'Sullivan ct al. 2003) and Chandra 
(Forman et al. 2002; Dosaj et al. 2002; Scharf et al. 2004; this work). The temperature of diffuse 
gas in NGC 1404 was found to be approximately isothermal with kT ~ 0.6 keV. The higher tem- 
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Fig. 1.— The 0.3 - 2.0 keV ACIS-S X-ray (left) and DSS (right) images of NGC1404 (SE, lower 
left) and NGC1399 (NW, upper right) galaxies in the Fornax cluster. The X-ray images have 
been binned by 2 (0".984 pixels), background subtracted, exposure corrected and coadded in WCS 
coordinates before smoothing with a 2" Gaussian. The image color scale has logarithmic stretch 
from 10~^ photons s~^ cm~^ arcsec"^ to 6 x 10~^ photons s~^ cm~^ arcsec"^ . The distance scale bar 
is 5 arcmin. The bright elliptical galaxy NGC 1404 exhibits a surface brightness discontinuity along 
its northwestern edge at a radius of ^ l'.4(8kpc) from its center. 

perature {kT ~ 0.75 keV) found by Loewenstein et al. (1994) was most likely a result of the large 
radius (3') aperture used, due to the ASCA point spread function, that extended well beyond the 
effective radius of NGC 1404. With the much better ROSAT PSPC resolution, Jones et al. (1997) 
found an asymmetric surface brightness distribution strongly suggesting that NGC 1404 was un- 
dergoing ram pressure stripping. Also using ROSAT PSPC data, Paolillo et al. (2002) measured 
the surface brightness asymmetry between the leading (northwest) sector and trailing (southeast) 
sector of NGC 1404, showing the steep gradient in the surface brightness in the northwest sector 
at r ~ 90" (8kpc); while O'Sullivan et al. (2003) observed a temperature jump for r > 2', across 
the leading edge of the galaxy and the surrounding Fornax ICM. 

In Figure 1 we present the 0.3—2 keV coadded image from three Chandra observations, totalling 
134.3 ks, (detailed in Section 2) of NGC 1404 falhng into the dominant cluster elhptical NGC 1399 
(left panel). For comparison we also show the Digitized Sky Survey optical image matched in WCS 
coordinates (right panel). NGC 1399 is located in the northwest (upper right) corner of each panel; 
while the bright elliptical NGC 1404 is located ~ 54kpc (9'. 8) in the plane of the sky to the southeast 
(lower left corner). The X-ray images have been individually background subtracted and corrected 
for telescope vignetting and spatial efficiency variations by means of exposure maps generated with 
standard CIAO tools assuming a fixed spectral energy of 0.9 keV. We clearly see the sharp edge in 
the surface brightness on the northwest side of NGC 1404 (in the direction of NGC 1399), formed by 
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gas bound to NGC 1404 shaped by the ram pressure of the cluster gas, and also the characteristic 
tail of ram-prcssurc-stripped gas to the galaxy's southeast. We use data from these three Chandra 
observations to measure the temperature and density of gas on either side of the surface brightness 
edge and calculate the dynamical motion of NGC 1404 relative to the Fornax ICM. This is an 
extension of earlier work on the Fornax cluster by Dosaj et al. (2002) using the Chandra ACIS-S 
observation (Obsid 319) to estimate the motion of NGC 1404 towards NGC 1399, and complements 
the recent, similar work by Scharf et al. (2004) using ACIS-I observations. Both assumed that the 
motion was primarily in the plane of the sky and, from the jump in gas pressure, obtained infall 
velocities of 466 ± SSkms"^ (Dosaj et al. 2002) and 660 ± 260 kms"^ (Scharf et al. 2004). 

Our discussion is organized as follows: In Section 2 we describe the observations and our 
data reduction and processing procedures. In Section 3 we discuss our analysis method and main 
results. In 3.1 we use spectral analyses to determine the temperature inside the surface brightness 
edge for NGC 1404 and outside the edge in the surrounding Fornax ICM showing that gas within 
the surface brightness edge is a factor ~ 3 cooler than the surrounding cluster gas. In Section 3.2 
we use coadded Chandra images from the three observations of NGC 1404 to extract the X-ray 
surface brightness profile across the edge and model the geometry of the cold front. We combine 
these data with the spectral information to construct the corresponding electron density profile 
across the edge, and in Section 3.3 we calculate the velocity of NGC 1404 relative to the Fornax 
ICM. In Section 4 we briefly summarize our results. Unless otherwise indicated all errors are 90% 
confidence levels and coordinates are J2000. Assuming a Hubble constant Hq = 75 km Mpc~^ in 
the currently preferred flat ACDM cosmology (17^ = 0.3, JIa = 0.7) and a redshift z = 0.004753, 
the redshift of the dominant elliptical galaxy NGC 1399, as representative for the Fornax cluster, 
the luminosity distance to the cluster is 19Mpc (Paolillo et al. 2002) and 1' corresponds to a 
distance scale of 5.49 kpc. 



2. Observations and Data Reduction 

Our data consist of three observations of galaxies NGC 1404 and NGC 1399 in Fornax taken 
with the Chandra X-ray Observatory using the Advanced CCD Imaging Spectrometer array (ACIS, 
Garmire et al. 1992; Bautz et al. 1998): a 57.4 ksec exposure taken on 18-19 January 2000 (when 
the focal plane temperature of the instrument was —110° C) using chips S2 and S3 (Obsid 319), a 
29.6 ks exposure (of NGC 1404 alone) taken on 13 February 2003 using the S3 chip (Obsid 2942), 
and a 47.3 ks exposure taken on 28-29 May 2003 using the ACIS-I array (Obsid 4174). For Obsid 
2942 and Obsid 4174 the focal plane temperature was —120° C throughout the observation. Each 
CCD chip is a 1024 x 1024 pixel array where each pixel subtends 0".492 x 0".492 on the sky. 

The data were analyzed using the standard X-ray processing packages, CIAO 3.0.2, FTOOLS 
and XSPEC 11.2. We corrected for the charge transfer inefficiency on the frontside-illuminated chips 
caused by exposure to low energy protons during a passage of the telescope through the Earth's ra- 
diation belts shortly after launch (Prigozhin et al. 2000), as well as for the time-dependent declining 



-5- 



efHciency of the ACIS detector due to the buildup of contaminants on the optical filter (Plucinsky 
et al. 2003) that is important at energies below 1.5 keV. Filtering removed events with bad grades 
(1, 5, and 7) and any events falling on hot pixels. Periods of anomalously high background (flares) 
were identified and were removed from the data along with periods of anomalously low count rates 
at the beginning and end of each observation. This resulted in useful exposure times of 56, 241 s, 
29,299 s, and 41,478 s, respectively, for the observations. 

Backgrounds for these observations were created from the source free ( blank sky) background 
sets appropriate for the CCDs, dates of observation and instrument configurations of interest 
(Markevitch et al. 2001)^ These were the merged 201,452 s (114,886 s) Period B source free 
dataset for the S2(S3) chip with ACIS-S at the aimpoint for Obsid 319, the merged 450,000 s 
Period D source free dataset for the S3 chip with ACIS-S at aimpoint for Obsid 2942, and the 
merged 550, 000 s Period D source free data set for the ACIS-I array with ACIS-I at the aimpoint 
for Obsid 4174. Identical cleaning, energy and spatial filters were applied to source and background 
data throughout, so that the normalization of the background is set, in each case, by the ratio of 
exposure times. We checked this normalization by comparing count rates between the source and 
background files in the 9 — 12keV range where particle background dominates. The source count 
rates were found to differ by < 3.5% from those of the blank sky sets for Obsid 319 and 2942, 
but were anomalously high (~ 16%) for Obsid 4174. We normalized the blank sky sets for that 
observation to agree with the observed 9 — 12keV source count rates. Thus we adopt 3.5% as 
the remaining uncertainty in the background levels for the three data sets. We identified point 
sources in the 16'. 8 x 16'. 8 field of the combined data in the 0.3 — lOkeV energy band using a 
multiscale wavelet decomposition algorithm set with detection threshold at 5a. The resulting 220 
source identifications were excluded from the following spectral and surface brightness analyses. 
We refer the reader to Angelini et al. (2001) for a detailed analysis of the point source population 
of this system. 



3. Results 

In Figure 2 we present surface brightness and temperature maps of the galaxy that clearly 
show the general structure of the X-ray emission from NGC 1404 and the surrounding ICM, i.e a 
sharp, bow-like shaped drop in surface brightness that is correlated to a similarly-shaped jump in 
temperature. In the northwest direction (the leading edge) the bow-shaped edge is nearly circular; 
while at larger angles the emission merges into a comet-like tail. In order to isolate the leading 
edge, we trace the surface brightness discontinuity with an elliptical segment centered at WCS 
coordinates (3''38"'52.1^ -35°35'39.6") with major(minor) axes of 8.3(8.0) kpc and position angle 
19.94°. The center of this edge-bounding ellipse is offset by < 2". 4 to the east of NGC 1404's optical 
center. We construct the surface brightness profile for the exposure-corrected, coadded image using 



^see http://cxc.harvard.edu/contrib/maxim/acisbg 



Fig. 2. — Surface brightness (left) and temperature map (right), matched in WCS coordinates, for 
NGC 1404 and the surrounding region. 

elhptical annuhi concentric to the edge-bounding ehipse and constrained to he within an angular 
sector extending from 35° to 91.4° measured counterclockwise from east (See Figure 3). For the 
surface brightness profile, radial intervals are chosen to increase(decrease) with logarithmic step 
size 1.05 as the radial distance is increased (decreased) from the position of the bounding ellipse. 
We need to obtain corresponding temperatures for gas in these regions. However, in order to obtain 
sufficient statistics for spectral analysis, we increase the logarithmic step size outside the bounding 
ellipse to 1.16. These spectral extraction regions are shown in Figure 3. We interpret region Wl as 
representative of the undisturbed, free streaming Fornax cluster gas and region WO as cluster gas 
directly outside the edge (Vikhlinin et al. 2001). In order to determine the spectral properties of 
gas inside the edge we use the elliptical annulus denoted by E in Figure 3 to maximize statistics 
in order to sharpen our determination of the metal abundance in NGC 1404 as well as the gas 
temperature inside the edge. 



Fig. 3. — Regions used in the edge analysis. The surface brightness profiles are extracted in elliptical 
annulli concentric to those shown and lieing within the angular sector shown (from 35° to 91.4° 
measured counterclockwise from East). Spectral extraction regions WO and Wl for the cluster 
gas outside the edge are also constrained to lie within that angular sector; while the full elliptical 
annulus E is used for the spectrum of NGC 1404 just inside the edge. 

3.1. Gas Temperatures 

Spectra from both Obsid 2942 and Obsid 4174 were fit simultaneously with XSPEC 11.2 
over the 0.3 — 5 keV energy range using APEC emission models for collisionally-ionized diffuse 
gas (Smith et al. 2001) corrected for absorption using Wisconsin photo-electric cross-sections 
(Morrison & McCammon 1983). Obsid 319 was not used in the spectral analysis due to calibration 
uncertainties for Period B observations. Counts were grouped with a pre-defined grouping resulting 
in channels of approximately constant logarithmic width. The hydrogen column density was initially 
fixed at its Galactic value 1.45 x lO^'' cm~^ (see http:/ /heasarc. gsfc.nasa.gov/, Archives & Software, 
nH:Column Density) with the gas temperature and abundance taken as free parameters. The 
hydrogen absorbing column was then allowed to vary to check the stability of the fit. The data 
showed no need for increased absorption either in the cluster gas or in the outer elliptical region E 
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Fig. 4. — (left) Best simultaneous two temperature APEC model fit to Obsid 2942 and Obsid 
4174 data for tfie elliptical annulus (E shown in Figure 3) just inside the edge where the second 
component models the cluster foreground emission, (right) Best simultaneous single temperature 
APEC model fit to Obsid 2942 and Obsid 4174 data for the Fornax cluster gas in the free stream 
region (Wl shown in Figure 3). Absorption is fixed at the Galactic value (nn = 1-45 x 10^° cm~^ ) 
in each region. 

in NGC 1404 just inside the edge. 

The results of the spectral analyses are summarized in Table 1. The best fits for the cluster free- 
stream region (Wl) and for the galaxy emission inside the edge (E) are shown in Figure 4. We find a 
temperature (1.531q'J3 keV) and abundance (0.42]']q^3 Zq) for gas in the free stream region. These 
values are in excellent agreement with XMM-Newton measurements by Buote (2002) for diffuse gas 
~ 8' from NGC 1399 (~ 2' from NGC 1404), who finds the data at large radius from NGC 1399 
well represented by a single temperature APEC model with temperature kT = 1.35 ± 0.15 keV 
(errors la) and abundance A 0.4 — 0.6 Zq. They arc also in reasonable agreement with Jones 
et al. (1997), who find kT = 1.26l|j:?^ keV and A = O.bltoHZQ for radfi r ~ 8' - 10' from 
NGC 1399 and by O'Sullivan et al. (2003), who find at the same radius kT ~ 1.2 - 1.45 keV and 
^ ~ 0.4 - 0.7^0 (see their Figure 1) using ROSAT PSPC data, and with Scharf et al. (2004) who 
find A;r ~ 1.6 using Chandra ACIS-I data. 

Due to limited statistics we are unable to freely vary all the parameters of the fit for region 

WO, the diffuse gas just outside the edge containing the stagnation region. We fix the abundance 
at the best fit value for region Wl and allow the temperature to vary. We find a temperature 
1.661q'^5 ^'^^ consistent with that for region Wl given the large errors. The somewhat cooler 
temperature (~ IkeV) measured by O'Sullivan et al. (2003) with ROSAT data when using a 
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circular annulus with mean radius 2' (llkpc) centered on NGC 1404 is most likely due to the 
fact that such a circular region contains cooler ram-pressure stripped gas from the trailing side of 
NGC 1404 (the debris tail) as well as the Fornax ICM. 

A single temperature APEC model provides a poor fit to the data inside the edge, as shown in 
Table 1. This is not surprising since a second component is warranted to model the contribution 
from foreground cluster emission. We model this emission with a separate APEC component with 
abundance and temperature fixed at the best fit cluster values. We find a temperature for the 
galactic diffuse gas inside the edge of 0.551q q2 keV with an abundance O.TSlg le ■^o- The fit is 
unchanged if the cluster component temperature is allowed to vary, although the errors on the 
fitted cluster temperature (1.48^0 22 keV) are large. The temperature we measure for diffuse gas in 
NGC 1404 is in good agreement with previous measurements given the differences in the models 
and spectral extraction regions. Using ROSAT PSPC data, Jones et al. (1997) found a mean 
temperature for NGC 1404 of 0.65to keV; while O'SuUivan et al. (2003) found a temperature 
profile consistent with isothermal with kT = 0.6 it 0.01 for radii r < 1'.6, similar to that found by 
Scharf et al. (2004) using ACIS-I data. We find, however, an abundance A significantly higher 
than that found by most previous authors {A ~ 0A4:Zq, Locwcnstein et al. 1994; A ~ 0.16 Zq, 
Jones et al. 1997; A ~ 0.35, O'Sullivan et al. 2003), and more consistent with solar or super- 
solar abundances expected for elliptical galaxies (Buote 2002; Brighenti &: Mathews 1999). This 
discrepancy is probably due to a combination of limited statistics, large heterogeneous extraction 
regions, and the use of simple single temperature models in those previous studies. We are in 



Table 1. APEC Model Spectral Fits 



Region 


RiruR-out 

kpCjkpc 


Source II(III) 
counts 


kT 
keV 


A 

Ze 


A 

lO-^^erg cm^s--*^ 


xVdof 


NGC 1404 E 

1 component 

2 component 

Fornax WO 
Fornax Wl 


5.1,8.2 

8.2,13 
13,21 


4071(2790) 

333(328) 
930(710) 


0.58 ±0.01 
1.53^ 

1 66+°-25 

1 C-q+0.10 

l.OO_o,i3 


0.32 ± 0.05 
0.42^ 

n 7Q+0.65 
"3-0.16 

0.42^ 
42+°-2 

^•^^-0.13 


1.0 
1.9 
0.846 


337/187 

276/185 
11.2/20 
104.5/69 



Note. — Column 2 lists the (inner, outer) radii of each extraction region shown in 3 measured from 
the center of the bounding ellipse. Source II (III) are the background subtracted source counts in 
the 0.3 — 5keV energy band for Obsid 2942 (Obsid 4174), respectively. Superscript / denotes a fixed 
parameter. A denotes the X-ray emissivity (cooling function) derived from XSPEC for the 0.3 — 2 keV 
energy band. Errors are 90% confidence limits. 
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agreement with the abundance measurement {A = O.TIq 3) for NGC 1404 found by Buote & Fabian 
(1998) using a multi-temperature MEKAL model to describe the ASCA data. 

In Table 1 we also list the X-ray emissivity (cooling function) A derived from the XSPEC best 
fit models for the energy range (0.3 — 2.0 kcV), the energy band we used for the surface brightness 
profiles in Section 3.2. In contrast to cluster analyses where the X-ray emissivity is dominated 
by the continuum free-free emission and depends only weakly on temperature and abundance 
(Vikhlinin et al. 2001), the emissivity of cool, metal-rich gas in galaxies is dominated by line 
emission and in this case varies by more than a factor of 2 between X-ray emitting gas in the 
cluster and that in the galaxy. Since the X-ray surface brightness depends on the product of the 
emission measure and the X-ray emissivity, high quality abundance measurements in galaxies are 
crucial for determining accurate electron densities from surface brightness profiles and subsequently 
constraining the dynamical motion of the system. 



3.2. Fitting the Surface Brightness Edge 

In Figure 5 we plot the X-ray surface brightness profile for NGC 1404, constrained by the 
wedge geometry shown in Figure 3 to sample only the region across the leading edge of the galaxy, 
as a function of the mean projected distance from the center of the edge-bounding ellipse. Three 
distinct regions are apparent. For small radii the data are well described by a power law with a 
break at 4 kpc. For radii 4 kpc ^ t < 8 kpc we see the surface brightness profile characteristic of a 
sharp edge, where the surface brightness decreases by an order of magnitude over a radial range of 
^ 1 — 2 kpc. This is characteristic for motion of a spheroid through a uniform gas and is consistent 
with a jump- like density discontinuity at the boundary. Outside the edge (r > 8 kpc) we see a 
surface brightness profile consistent with cluster gas with a possible enhancement at r ~ 10 kpc, 
near the edge. The gradual rise in surface brightness for r > 20 kpc (< 34 kpc from NGC 1399) is 
due to halo emission from NGC 1399 (See Figure 6). 

In order to fit the surface brightness in Figure 5, we need to model the form of the electron 
density in the galaxy NGC 1404 and the Fornax ICM. We assume, as in Vikhlinin et al. (2001), 
that the cold front inside the edge can be described by an isothermal sphere whose density is well 
modelled by a power law 

(f \ —a 
j (1) 
''edge ' 

where Tcdgc is the position of the edge and Agdge is the electron density in the cold front inside the 
galaxy at the position of the edge (r = redge)- 

The density model for the cluster emission is more problematical since it requires knowledge of 
the density profile across the full extent of the cluster along the line of sight (r < 1.5 Mpc) and may 
contain local variations due to the presence of significant STibclustcr mass concentrations common 
in dynamically young clusters (see, e.g. Schindler et al. 1999 for Virgo; Machacek et al. 2002 
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Fig. 5. — (left) Surface brightness profile across the edge from elliptical annulli bounded by the 
sector shown in Figure 3 as a function of the distance from the center of the bounding ellipse. The 
models shown represent power law density distributions inside the edge (4kpc < r < redge) and 
two beta model fits to the Fornax ICM (r > Tcdgc): P = 0.35, = 45" from Jones et al. (1995) 
(PL+ICMl, solid hnc) and /? = 0.3, = 394" (PL+ICM2, dashed line) from a three component 
beta model fit to NGC 1399. The dotted line shows the best fit for a single power law (SPL) 
model for radii 1 kpc < r < redge inside the edge and the ICM2 model outside the edge, (right) 
Corresponding electron density profiles for PL+ICMl (solid line) and PL+ICM2 (dashed line) 
models showing the electron density discontinuity across the edge. 



and AbdelSalam et al. 1998 for Abell 2218; Drinkwater et al. 2001 for Fornax). We adopt simple 
isothermal beta models for the surface brightness and gas density of the Fornax cluster gas outside 
the edge such that the surface brightness S{r) is given by 

/ / r \ 2n -3/3+0.5 

S(r) = So (! + (-)) . (2) 

where S{r) is the surface brightness at a given radius, is the central surface brightness, and Tc 
is the core radius. Then the corresponding electron density is given by 

ne(r)=no(l+(^)')~'^^'. (3) 

In Figure 6 we show two multiple-component beta model fits to the surface brightness profile 
of NGC 1399, with point sources and emission from NGC 1404 excluded, where the fit to the data 
is the sum of individual components given by Equation 2. In each case the innermost central region 
is modelled with Vc = 6" and (5 = 0.65. This result is in good agreement with the fit to NGC 1399 
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Fig. 6. — Beta model fits to the surface brightness profile of NGC 1399 as a function of the radial 
distance from the center of NGC 1399: (left) Two component fit with (rc,/?) given by (6", 0.65, 
dotted line) and (45", 0.35, dashed line) for the core and outer regions, respectively, (right) Three 

component fit with Vc, (3 for the core, halo, and cluster given by (6", 0.65, dotted line), (98", 0.5, 
dashed hne), and (394", 0.3,long dashed line), respectively. The solid line denotes the sum of the 
components in each case. 

by O'Sullivan et al. (2003), who find rc = 6.6 d= 0.6" and P = 0.59 ± 0.1, although they use an 
elliptical beta model with axis ratio 1.23 rather than our spherically symmetric one. 

Surface brightness profiles at larger radii have been measured by Jones et al. (1995) and 
Paolillo et al. (2003) using ROSAT PSPC data. Jones et al. (1997) found the surface brightness 
profile beyond the core of NGC 1399 well represented by a single beta model with core radius 
rc = 45" and (3 = 0.35 out to a radius of 40'(220kpc). Paolillo et al. (2003) modelled the surface 
brightness profile of NGC 1399 and the Fornax ICM with a three component beta model, whose 
third (cluster) component dominated emission beyond 7'. Their model is consistent with the model 
of Jones et al. (1997) for radii r < 200 kpc, where the ROSAT PSPC corrections are well known, 
but falls off more rapidly at larger radii. 

Similarly, we find our data equally well described by either a two or three component beta 
model over the radial range of our observation. The left panel of Figure 6 models the data with a 
two component beta model where the outer component is parameterized by rc = 45" and P = 0.35, 
hereafter called ICMl, as found by Jones et al. (1997). The right panel shows a three component 
beta model fit to the same data, where the surface brightness at large radii is described by the 
sum of a central component with (vc, f3) the same as above and two additional components, the 
galactic halo gas with (rc = 98", (3 = 0.5) and cluster gas with (rc = 394", /3 = 0.3), hereafter called 
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ICM2. The surface brightness edge in NGC 1404 is located at a projected distance of ^ 46kpc 
from the center of NGC 1399 where the second (third) beta model component in the two (three) 
component fit dominates. Since this represents the minimum physical distance between the edge 
and the center of the dominant cluster elliptical NGC 1399 (and the overlap between components 
in this region is small) , we adopt the density beta model from Equation 3 with parameters fixed by 
those components as representative of the Fornax cluster gas in our analysis of the edge. Our data 
do not probe cluster distances large enough to distinguish between these two density models. Thus 
we consider both models in order to investigate how uncertainties in the cluster density distribution 
at large distances affect our results. 

In order to extract the density distribution from the surface brightness profiles, we need to 
make additional assumptions about the relative three dimensional geometry of NGC 1404 with 
respect to NGC 1399. The comet-like tail to the southeast of NGC 1404 seen in Figures 1 and 
3 is evidence for ram pressure stripping of gas from NGC 1404. Since ram pressure stripping is 
more likely in the denser, central regions of the cluster (i.e. near the dominant elliptical galaxy 
NGC 1399), we assume that NGC 1404 and NGC 1399 both lie near the plane of the sky such that 
the projected distance 9'. 8 is representative of their physical radial separation. The direction of 
the tail, pointing to the southeast and away from the direction of motion, indicates that NGC 1404 
is moving to the northwest toward NGC 1399. Finally, the sharpness of the edge suggests that 
we are able to directly view the stagnation region, such that the inclination of the motion with 
respect to the plane of the sky should not be too large (Mazzotta et al. 2001). These are the same 
assumptions used in the analysis of the infalling subcluster in Abell 3667 (Vikhlinin et al. 2001). 
Our analysis method closely follows that work. 

We calculate the surface brightness by integrating the square of the density times the cor- 
responding X-ray emissivity, given in the sixth column of Table 1, along the line of sight. The 
data are fit using a multi-variate chi-square minimization scheme allowing the position of the edge 
(redge)) density power law index inside the edge (a), and galaxy electron density at the edge (^edge) 
to vary. We also allow the cluster density normalization (no in Equation 3) to be a free parameter 
to partially compensate for possible local density fluctuations in the cluster ICM when extracting 
the cluster electron density near the edge. 

In Figure 5 we show the results for three density models: a power law (PL) density model 
for the edge region (4 kpc < r < redge) interior to NGC 1404 with the cluster density modelled 
by either ICMl or ICM2, and a single power law (SPL) model for r > 1 kpc, with no break at 
r ~ 4 kpc. The ICM2 model is used for the cluster density in this latter case. As the dotted line 
in Figure 5 shows, the surface brightness profile for radii < 4 kpc is well reproduced by this simple 
power law density model for the galaxy (SPL-I-ICM2) with a best fit power law index a = 1.54 
(corresponding to a beta model with /? = 0.47 for r >> Tc), in reasonable agreement with the beta 
model results (/? = 0.51, Tc = 0.57 kpc) for NGC 1404 by Paolillo et al. (2002) and (/3 ~ 0.5) by 
Scharf et al. (2004). However, this model is unable to reproduce the shape of the surface brightness 
distribution in the edge region (4 kpc < r < redge)- Both density models PL-MCMl (solid line) 
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and PL+ICM2 (dashed line) reproduce the surface brightness profile across the edge well. The fit 
parameters for these models are summarized in Table 2. The electron density we derive by fitting 
the surface brightness profile directly for gas inside the edge (3.9 x 10~'^ cm^^ for the PL+ICMl 
model and 4.3 x 10"'^ cm"'^ for the PL+ICM2 model) are in excellent agreement with the beta 
model result (~ 4 x 10~^ cm~^ ) result by Scharf et al. (2004) , although our value for the cluster 
gas is ~ 20% higher. 

3.3. Bounding the Galaxy's Dynamical Motion 

Following Vikhlinin et al. (2001) we use the densities and the temperatures from our surface 
brightness and spectral fits to calculate the gas pressures in the free stream region and in the galaxy 
gas just inside the edge, which is in pressure equilibrium with gas at the stagnation point (where 
the relative cluster gas velocity vanishes). The difference between the two pressures is interpreted 
as the ram pressure of the intracluster medium on the leading edge of the galaxy (the cold front); 
while the ratio of the pressures allows us to determine the Mach number Mi = v/ci (where ci is 
the speed of sound in the cluster free stream region) for the cold gas cloud moving through the 
hot ICM. For completeness, the pressure-jump versus Mach number relation, plotted in Figure 7, 
is given by : 



Table 2. Best Fit Edge Position and Electron Densities 



Property 


PL+ICMl Model 


PL+ICM2 Model 


redge (kpc) 


8.01 


7.64 


a 


0.84 


0.54 


Win (10~^ cm~^ ) 


3.9 


4.3 


nont (10-3 cm-3) 


0.82 


0.72 



Note. — Tedge is the position of the edge measured from the 
center of the bounding ellipse in Figure 3, a is the power law 
index of the spherically symmetric electron density distribu- 
tion inside the edge, and nin (jT-out) is the electron density just 
inside (outside) the edge. 
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Mach number 

Fig. 7. — Pressure jump as a function of Mach number as given by Equation 4 for Mi < 1 and 
Equation 5 for Mi > 1 (dot-dashed hne). The horizontal hues denote the value of the pressure 
jump for the PL+ICMl model (solid line) and for the PL+ ICM2 model (dashed line) fits to the 
surface brightness profile across the edge. 



for Ml < 1 (subsonic regime); while 
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for Ml > 1 (supersonic regime), with 7 = 5/3 the adiabatic index for a monatomic, ideal gas 
(Landau &: Lifshitz, 1959). 

In Table 3 we summarize the results of this analysis. From Figure 7 we find a Mach num- 
ber for the galaxy gas cloud moving through the cluster gas of 0.83 (1.03) for the PL+ICMl 
(PL+ICM2) best fit density models, respectively. The sound speed for the completely ionized 
(mean molecular weight // = 0.6) intracluster gas at a temperature 1.531q{3 stream 
region is ci = 638kms~^. Thus the speed v = Mici of NGC 1404 relative to the ICM is 
found to be u = 531(657) kms for these two models, in excellent agreement with the estimate 
(v ~ 660 lb 260 km s^-*^ ) by Scharf et al. (2004). If we assume that the Fornax cluster ICM, at least 
that part of it that impacts the leading edge of NGC 1404, is at rest relative to NGC 1399 and use 
the radial (line of sight) velocities 1947 km s~^ for NGC 1404 and 1425 km s"^ for NGC 1399 given 
in NED (Graham et al. 1998), then the radial velocity between these galaxies would be 522 kms~^ . 
The relative transverse component of the velocity between NGC 1404 and the Fornax ICM is then 
98 (399) kms^^ , respectively, implying an inclination angle for the motion with respect to the plane 
of the sky of ^ ~ 79° (53°), respectively, for our two models. However, Mazzotta et al. (2001) stress 
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that the sharpness of the edge is a strong function of the inchnation angle ^. Given the sharp- 
ness of the observed edge in NGC 1404, these derived inchnation angles arc uncomfortably large. 
If we adopt instead the optical relative radial velocity Vr = 450 km s^"*^ measured by Drinkwater 
et al. (2001), who identify NGC 1404 as part of a larger, high velocity, infalling Fornax subcluster, 
the inferred transverse velocities for NGC 1404 increase to ~ 292 (479) km s^^ , respectively, for 
our two models. While improved, the inferred inclination angle of the motion with respect to the 
plane of the sky, ^ = 58° (43°), is still large. 

There are several possible solutions to this dilemma. First, NGC 1404 and NGC 1399 may 
not both lie in the plane of the sky. NGC 1404 may be undergoing a fiy-by of the dominant 
elliptical, with a non-negligible impact parameter rather than a direct collision, such that the use 
of the projected distance between the galaxies underestimates the physical radial distance between 
them and overestimates the value of the cluster electron density at the edge. Thus the density 
jump should be viewed as a lower bound on the physical density discontinuity, Mach number, and 
relative velocity of NGC 1404 with respect to the intracluster gas. A second possibility is that 

Table 3. Derived Velocities and Edge Properties for Two Cluster ICM Models 



Property 


PL+ICMl Model 


PL+ICM2 Model 


neo/nei 


4.78 


6.04 


To/Ti 


0.36 


0.36 


Po/pi 


1.7 


2.1 


Ml 


0.83 


1.03 


v(kms~^) 


531 


657 


ft(kms~^) 


98(298) 


399(478) 


^ (degrees) 


79(58) 


53(43) 



Note. — ng, T, and p denote the electron density, 
temperature and pressure in each region specified by 
subscripts and 1 for galaxy gas just inside the edge and 
cluster gas in the free stream region, respectively (See 
Vikhlinin et al. 2001). Mi is the Mach number, v is the 
relative total velocity, and Vt and ^ are the component of 
velocity in the plane of the sky and inclination angle of 
motion with respect to the plane of the sky, respectively, 
assuming a relative radial velocity between NGC 1404 
and NGC 1399 of Vr = 522 km s"^ from NED {vr = 
450kms~^ from Drinkwater et al. 2001). 
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NGC 1399 is moving with respect to the Fornax ICM and may, itself, be experiencing ram pressure 
stripping due to this motion (PaoHllo et al. 2002). Then the optically measured relative radial 
velocity between the two galaxies may not accurately reflect the relative radial velocity between 
NGC 1404 and the surrounding cluster gas. 

Despite these uncertainties it is useful to consider what our analysis implies for the increased 
temperature and density expected in the stagnation region. For the subsonic motion (M = 0.83) 
found for model PL+ICMl and using the cluster free stream as the reference region, Bernoulli's 
equation predicts a density increase of 37% leading to a surface brightness enhancement by a 
factor of ~ 1.9. This is consistent with the possible ^ 70% surface brightness enhancement seen 
justoutside the edge in Figure 5. The tempcratiue in the stagnation region would be expected to 
rise by ~ 23% by adiabatic compression to ~ 1.88 keV. However, the stagnation region is small 
(of the order of a galaxy diameter) such that the observed temperature would be lower because 
of dilution by cooler cluster gas along the line of sight. Our measurements for the temperature 
in this region hint of a possible temperature rise consistent with such adiabatic heating, but our 
errors are too large because of limited statistics for a definitive measurement. For model PL+ICM2 
(M 1.03) the temperature and density behind the shock would be expected to increase by ~ 3% 
and ~ 5%, respectively, with temperature and density in the stagnation region similar to that for 
model PL+ICMl. Such a weak shock would not be expected to be seen in our data. 



4. Conclusions 

In this paper we have used data from three Chandra observations using the ACIS array (a 
57.4 ksec exposure with chips S2 and S3 taken on 18-19 January 2000, a 29.6 ks exposure using chip 
S3 taken on 13 February 2003 , and a 47.3 ks exposure using the ACIS-I array taken on 28-29 May 
2003) of the elliptical galaxy NGC 1404 falling toward the dominant elliptical galaxy NGC 1399 in 
the Fornax cluster to bound the relative physical motion of NGC 1404 with respect to the Fornax 
cluster ICM. We find: 

1. The temperature and abundance of the galaxy gas just inside the edge in NGC 1404 is 
0.55lo 02 ^^'^ O-T^^lo 16-^0) typical of elliptical galaxies. Gas temperatures for the cluster 
gas in the free-stream region (Wl) and region next to the edge containing the stagnation point 
(WO) are 1.53l|j:^^ keV and 1.66l^;^| keV, respectively, with a cluster abundance 0.42j:°;^3ZQ 
measured in Wl. 

2. The X-ray emissivity in NGC 1404 is more than a factor 2 greater than that of the surrounding 
cluster gas. Since the X-ray emissivity is dominated by line emission in cool, metal rich 
galaxies, high quality abundance measurements in galaxies are crucial for determining the 
density discontinuity across the edge and constraints on the dynamical motion of the galaxy 
through the ICM from these edge analyses. 
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3. The surface brightness profile in NGC 1404 in the direction towards NGC 1399 obeys a power 
law for 1 kpc < r < 4 kpc, a break at ~ 4 kpc, and a sharp edge where the surface brightness 
decreases by more than an order of magnitude over 1 — 2 kpc. Outside the edge the surface 
brightness profile is consistent with cluster gas with a possible 50 — 70% enhancement at 
r ~ 10 kpc. 

4. The best fit edge position is found to be at 7.6 — 8.0 kpc from the center of the edge-bounding 
ellipse with a weak dependence on the large distance behavior of the cluster gas density. A 
single power law density model within NGC 1404 provides a poor fit to the data for r > 4 kpc. 
The surface brightness profile across the edge is well modelled with a power law electron 
density in NGC 1404 with index a = 1.54 (beta model index f3 = 0.47) for 1 kpc < r < 4 kpc, 
a break at 4 kpc after which a flattens to ~ 0.5 — 0.8 with a weak dependence on the 
beta model parameters assumed for the cluster gas. The sharpness of the edge is consistent 
with a jump- like density discontinuity with electron density in NGC 1404 inside the edge of 
3.9 — 4.3 X 10~^ cm~^ and in the cluster free stream region of 7 — 8 x 10~^ cm~^ . 

5. The pressure jump, po/pi ~ 1.7 — 2.1, between the cluster free stream region and gas at the 
edge implies a Mach number Mi ~ 0.83 — 1.03 relative to gas in the cluster free stream region 
and thus a velocity v ~ 531 - 657 km s~^ for NGC 1404 through the surrounding ICM. This 
motion predicts an enhancement in the surface brightness of the cluster gas in the stagnation 
region by a factor ~ 1.9 due to adiabatic compression, consistent with that seen. 

6. Using values for the relative radial velocity between NGC 1404 and NGC 1399 drawn from 
the literature as representative of the relative radial velocity between NGC 1404 and the 
surrounding cluster gas, the inferred value for the inclination angle of the motion with respect 
to the plane of the sky is large, > 40°, in apparent contradiction with the observed sharpness of 
the edge. This may signal the presence of a significant impact parameter between NGC 1404 
and NGC 1399 or that NGC 1399 is also moving radially with respect to the cluster gas. 
Thus, more detailed hydrodynamic modelling of this system is required to resolve the full 
three dimensional motion of this system. 
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